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ABSTRACT 



Precise radial-velocity measurements with the HARPS spectrograph reveal the presence of two planets orbiting the solar-type star 
HD 45364. The companion masses are m sin; = 0.187 M Jup and 0.658 M Jup , with semi-major axes of a = 0.681 AU and 0.897 AU, and 
eccentricities of e = 0.168 and 0.097, respectively. A dynamical analysis of the system further shows a 3:2 mean motion resonance 
between the two planets, which prevents close encounters and ensures the stability of the system over 5 Gyr. This is the first time that 
such a resonant configuration has been observed for extra-solar planets, although there is an analogue in our Solar System formed by 
Neptune and Pluto. This singular planetary system may provide important constraints on planetary formation and migration scenarios. 

Key words, stars: individual: HD 45364 - stars: planetary systems - techniques: radial velocities - methods: observational 



1. Introduction 

At present, about 25% of the known exoplanets are in multi- 
planetary systems. The older the radial-velocity planet-search 
surveys, the higher the fraction of multi-planet families detected 
in these programs. Taking into account the still strong bias 
against long-period and/or low-mass planet detection, and the 
enhanced difficulty of fully characterizing systems with more 
than one planet, the fraction of known multi-planet systems is 
certainly still a lower limit. It appears likely that a high number, 
if not the majority, of planet-host stars form systems of planets 
rather than isolated, single planetary companions. 

Among the known multi-planet systems, a significant frac- 
tion are in mean motion resonances, the majority of which are 
in low-order resonances. The 2:1 resonance is the most com- 
mon (HD 73526, HD 82943, HD 1283 1 1, GJ 876), but other con- 
figurations are observed as well, such as the 3:1 resonance in 
HD 75732 or the 5:1 resonance in HD 202206. These resonances 
most probably arise from evolutionary processes as migrating 
planets forming in the protoplanetary disc become trapped. In 
our Solar System, we also find mean motion resonances in the 
satellites of the giant planets, or between these planets and sev- 
eral asteroids. The most well-known example is the Io-Europa- 
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* Based on observations made with the HARPS instrument 
on the ESO 3.6 m telescope at La Silla Observatory un- 
der the GTO programme ID 072.C-0488. The table with the 
radial velocities is available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/????| 



Ganymede system in a 4:2:1 resonance, or the Neptune-Pluto 
system in a 3:2 resonance. While the satellites are believed to 
achieve resonant configurations after tidal evolution of their or- 
bits, the resonances between planets and asteroids were probably 
formed after the inward or outward migration of the planets dur- 
ing the early stages of the evolution of the Solar System. 

The presence of two or more interacting planets in a system 
dramatically increases our potential ability to constrain and un- 
derstand the processes of planetary formation and evolution. The 
dynamical analysis of such systems is then very useful, first for 
constraining the system evolution history and second for deter- 
mining the system "structure" in terms of orbital content. 

Multi-planet systems are naturally found in planet-search 
programs. However, improving the precision of the radial- 
velocity measurements greatly helps their detection and com- 
plete characterization. The HARPS search for southern extra- 
solar planet is an extensive radial-velocity survey of some 2000 
stars of the solar neighborhood conducted with the HARPS spec- 
trograph on the ESO 3.6-m telescope at La Silla (Chile) in the 
framework of the Guaranteed Time Observations granted to the 
HARPS building consortium dMavor et al.1 120031) . About half 
of the HARPS GTO time is dedicated to very high-precision 
measurements of non-a ctive stars select ed from the CORALIE 
planet-search program (lUdrv et al.ll2.000h and stars with already 
known giant planets, while searching for lower mass compan- 
ions. This prog ram reveals itself as very efficient in finding 
multi -Neptune (Lovis et al.l 12006) and multi-Super Earth sys- 
tems dMavor et alj|2008l) . 
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Table 1. Observed and inferred stellar parameters of HD 45364. 



Parameter 




HD 40307 


Sp 




KOV 


V 


[mag] 


8.08 


B-V 


[mag] 


0.72 


IT 


[mas] 


30.69 ± 0.81 


My 


[mag] 


5.51 


Tss 


[K] 


5434 ± 20 


logg 


Legs] 


4.38 ± 0.03 


[Fe/H] 


[dex] 


-0.17 ±0.01 


L 


[L B ] 


0.57 


M, 


[Ma] 


0.82 ± 0.05 


v sin i 


[km s" 1 ] 


1 


logtf HK 




-4.94 


ProtClogl^) 


[day] 


32 



Table 2. Orbital parameters for the two bodies orbiting 
HD 45364, obtained with a 3-body Newtonian fit to observa- 
tional data (Fig.[T]>. 



Photometric and astrometric data are from the Hip parcos catalogu e 
dESAll 997) and the stellar physical quantities from lSousa et 



From the HARPS high-precision survey, we present an in- 
teresting system of two planets in a 3:2 mean motion reso- 
nance around HD 45364, a configuration not observed previ- 
ously among extra-solar planet, but similar to the Neptune-Pluto 
pair in the Solar System. Section 2 gathers useful stellar infor- 
mation about HD 45364, its derived orbital solution is described 
in Sect. 3, and the dynamical analysis of the system is discussed 
in Sect. 4. Finally, conclusions are drawn in Sect. 5. 

2. Stellar characteristics of HD 45364 

The basic photometric (KOV, V = 8.08, B-V- 0.72) and as- 
trometric (n- 30.69 mas) pro perties of H D 45364 were taken 
from the Hipparcos catalogue dESAll 1997b . They are recalled in 
Table[T]with inferred quantities such as the absolute magnitude 
(My =5.51) and the stell ar physical characte ristics derived from 
the HARPS spectra by ISousa et all J2008). For the complete 
high-precision HARPS sample (including HD 45364), these au- 
thors provided homogeneous estimates of effective temper- 
ature (r eff = 5434 + 20 K), metallicity ([Fe/H] = -0.17 + 0.01), 
and surface gravity (log g — 4.38 + 0.03) of the stars. A projected 
low rotational velocity of the star, vsinz = 1 kms~' was derived 
from a calibration of the width of t he cross-correlatio n function 
used in the radial-velocity estimate dSantos et al.ll2002l) . 

HD 45364 is a non-active star in our sample with an activ- 
ity indicator log R' HK of -4.94. No significant radial-velocity 
jitter is thus expected for the star. From the activity indicator, 
we also derive a s tellar rotation period P ro , = 32 day (following 
iNoveset al.lll984l) . 

HD 45364 has a subsolar metallicity with [Fe/H ] = -0.17 un- 
like m ost of the gaseous giant-planet host stars dSantos et al.l 
|2004|) . According to simulations of planet formation based 
on the core-accretion paradigm, moderate metal-deficiency 
does not, however, pr event planet formation (llda & Lin 2004; 
Mordasini et al. 2007). Taking into account its subsolar metal- 
licitv. lSousa et al.l l2008) derived a mass of 0.82 M Q for the star. 
From the color index, the derived effective temperature, and the 
corresponding bolometric correction, we estimated the star lu- 
minosity to be 0.57 L . 



3. Orbital solution for the HD 45364 system 

The HARPS observations of HD 45364 started in December 2003 
and have now been going for about four years and a half. From 





[unit] 


HD 45364 b 


HD 45364 c 


Date 


[JD-2400000] 


53500.00 (fixed) 


V 


[km/s] 


16.4665 ± 0.0002 


P 


[day] 


226.93 ± 0.37 


342.85 ± 0.28 


A 


Tdesl 


105.76 ± 1.41 


269.52 + 0.58 






0.1684 + 0.0190 0.0974 + 0.012 




[deel 


162.58 ± 6.34 


7.41 ± 4.30 


K 


[m/s] 


7.22 + 0.14 


21.92 + 0.43 


i 


[cleg] 


90 (fixed) 


90 (fixed) 


a i sin; 


[10- 3 AU] 


0.1485 


0.6874 


f(M) 


[10- 9 M,,] 


0.0085 


0.3687 


M sin i 


[M Jup ] 


0.1872 


0.6579 


a 


[AU] 


0.6813 


0.8972 


N 




58 




Span 


[day] 


1583 




rms 


[m/s] 


1.417 








2.789 





Errors are given by the standard deviation a and A is the mean longitude 
of the date (A = a> + M). The orbits are assumed co-planar. 
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Fig. 1. HARPS radial velocities for HD 45364, superimposed on 
a 3-body Newtonian orbital solution (Tabled. 



the first stages of the observations, peculiar variations in the 
radial velocities (Fig.[TJ have shown the presence of one or 
more companions in the system. After 58 measurements, we 
are now able to determine the nature of these bodies. Using 
a genetic algorith m combined with the iterative Levenberg- 
Marquardt method dPress et al.lll992h . we first attempted to fit 
the complete set of radial velocities using a model with two 
Keplerian orbits. This fit yields an adjustment of yfx^ = 2.71 
and rms = 1.38 ms _I with two planets, one at P = 225.8 day, 
e = 0.174, and a minimum mass of 0.186 Mj up , and another at 
P = 343.9 day, e = 0.097, and a minimum mass of 0.659 M Jup . 

Due to the proximity of the two planets and their high min- 
imum masses, the gravitational interactions between these two 
bodies are strong. This prompts us to fit the observational data 
using a 3-body Newtonian model, assuming co-planar motion 
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Fig. 2. Radial velocity differences between the two independent 
Keplerian model and the 3-body Newtonian model (Tabled. 
Data coincides at JD=2453776 (Feb. 9 th 2006). We also plot 
the velocity residuals of the 3-body fit (Fig.Q]). With the cur- 
rent HARPS 's precision for this star, we expect to observe these 
differences within ten years. 



perpendicular to the plane of the s ky, similarly to wha t has been 
done for the system HD 202206 dCorreia et alj|2005l) . The or- 
bital parameters corresponding to the the best fitted solution are 
listed in Table|2] We get identical results for \Jx* and veloc- 
ity residuals as obtained with the two-Keplerian fit. Although 
there is no significant improvement in the fit, an important dif- 
ference exists: the new orbital parameters for both planets show 
some deviations from the two-Keplerian case. In particular, the 
arguments of the periastrons (at) show a difference of some de- 
grees. We then conclude that, despite being unable to detect the 
planet-planet interaction in the present data, the orbits undergo 
important perturbations, and we expect to detect this gravita- 
tional interaction in the future. In Fig.|2j we plot the two best 
fit models evolving in time and clearly observe detectable devi- 
ations between the two curves that appear within ten years. The 
3-body Newtonian fit provides a more accurate approximation of 
the HD 45364 planetary system, and the orbital parameters thus 
determined will be adopted as reference henceforward (Table|2]). 

We also fitted the data with a 3-body Newtonian model for 
which the inclination of the orbital planes, as well as the node of 
the outer planet orbit, were free to vary. We were able to find a 
wide variety of configurations, some with low inclination values 
for one or both planets, that slightly improved our fit to a mini- 
mum = 2.32 and rms = 1.14m/s. However, all of these de- 
terminations are uncertain, and since we also increase the num- 
ber of free parameters by three, we cannot say that there has 
been an improvement with respect to the solution presented in 
Table|2] The inclination of the orbits therefore remain unknown, 
as do their true masses. 

The residuals around the best fit solution are small, but re- 
main slightly larger than the internal errors (Fig.[T|i. We may 
then ask if there are other companions with different orbital pe- 
riods. For this purpose, we used a genetic algorithm, since we 
were unable to clearly isolate any significant peak in the fre- 
quency analysis of the residuals. The inclusion of additional 
companions in the system allows us to reduce the value of y/x* 
slightly, although this can be justified as a natural consequence 
of increasing the number of free parameters. Identical adjust- 
ments can be obtained with many orbital periods, as different 
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Fig. 3. Variation in the resonant argument, 8/, = 2Aj - 3A C + 
a>b (a) and in the secular argument, Aa> = a>b - a> c (b), with 
time. 6b is in libration around 0°, with a libration period Pi g st 
18.16yr, and a principal amplitude of about 68.4° (Table|4|. Aa> 
is in libration around 180°, with a libration frequency g^ = gi~ 
g2 (corresponding to a period P^ w =2 413.9 yr), and a maximum 
amplitude of about 36.4°. 



as 5 or 18 days, frequently with very high eccentricity values. 
Therefore, no other companion can be conclusively detected in 
the residuals from the orbital solution listed in Table|2] The best 
fit solution was obtained by adding a linear drift to the data, with 
slope = -0.86 + 0.09 ms"'/yr, allowing us to reduce the value 
of \Jx* to 2.42 and the rms = 1.21 m/s, while the orbital param- 
eters of the two planets remain nearly the same. The solution 
in Table|2]must then be considered to be the best determination 
achievable so far, and a longer tracking of the system will pro- 
vide more accurate orbital parameters for the HD 45364 system. 

4. Dynamical analysis 

We now briefly analyze the dynamics and stability of the plane- 
tary system given in Table[2] Due to the two planets' proximity 
and high values of the masses, we expect that both planets are 
affected by strong planetary perturbations from each other. The 
present orbits of the two planets almost cross (Fig.|6]l, and unless 
a resonant mechanism is present to avoid close encounters, the 
system cannot be stable. 

4.1. The 3:2 mean motion resonance 

The ratio between the orbital periods of the two planets deter- 
mined by the fitting process (Table|2| is P c /Pb = 1.511, sug- 
gesting that the system may be trapped in a 3:2 mean motion 
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0.62 0.64 0.66 0.68 0.70 0.72 0.74 0.84 0.86 0.88 0.90 0.92 0.94 0.96 

a b {AU) a c {AU) 

Fig. 4. Stability analysis of the nominal fit (Table© of the HD 45364 planetary system. For a fixed initial condition of the outer 
(a) and inner planet (b), the phase space of the system is explored by varying the semi-major axis and eccentricity of the 
other planet, respectively. The step size is 0.001 AU in semi-major axis and 0.005 in eccentricity. For each initial condition, the full 
system i s inte grated numeric ally over 10 kyr and a stability criterion is derived with the frequency analysis of the mean longitude 
dLaskarll 1 990l 1 1 993h . As in lCorreia et alj d2005l) . the chaotic diffusion is measured by the variation in the frequencies. The "red" 
zone corresponds to highly unstable orbits, while the "dark blue" region can be assumed to be stable on a billion-years timescale. 
The contour curves indicate the value of x 1 obtained for each choice of parameters. It is remarkable that in the present fit, there is 
perfect correspondence between the zone of minimal^ 2 and the 3:2 stable resonant zone, in "dark blue". 



Table 3. Fundamental frequencies for the orbital solution in 
Tabled 





Frequency 


Period 




°/yr 


yr 


rib 


576.429624 


0.624534 


n c 


384.033313 


0.937419 


gi 


-0.759309 


474.115 


g2 


0.110472 


3258.74 


k 


19.820696 


18.1628 



n b and n c are the mean motions, g\ and g2 are the secular frequencies 
of the periastrons, and l e is the libration frequency of the resonant angle 
6 b = 2A b - 3A C + oj b . Indeed, we have 2n b - 3n c + gj = 0. 

resonance. To test the accurancy of this scenario, we performed 
a frequency analysis of the orbital solution listed in Table©com- 
puted over 1 00 kyr. The orbits of the planets are integrated wit h 
the symplectic integrator SABAC4 of lLaskar & Robutell d2001h . 
using a step size of 0.02 years. We conclude that in the nomi- 
nal solution of Table© the two planets in the HD 45364 system 
indeed show a 3:2 mean motion resonance, with resonant argu- 
ment: 

b = 2A b - 3A C + co b . (1) 

The fundamental frequencies of the systems are the two mean 
motions n b and n c , the two secular frequencies of the periastrons 
gi and g2, and the libration frequency of the resonant argument 
le (Table©. These frequencies are not independent because, due 
to the 3:2 resonance, we have up to the precision of the determi- 
nation of the frequencies (w 10~ 10 ), 

2n b - 3n c + g\ = . (2) 

The resonant argument 6 b is in libration around 0°, with a 
libration period 2n/lg = 18.16 yr, and an associated amplitude 



of about 68.44 degrees (Fig.©, Table©. For the complete so- 
lution, the libration amplitude can reach more than 80 degrees 
because additional periodic terms are present. In Table© we pro- 
vide a quasi-periodic decomposition of the resonant angle 6 b in 
terms of decreasing amplitude. All the quasi-periodic terms are 
easily identified as integer combinations of the fundamental fre- 
quencies (Table©. Since the resonant angle is modulated by a 
relatively short period of about 18 years, the observation of the 
system over a few additional decades may provide an estimate 
of the libration amplitude and thus a strong constraint on the pa- 
rameters of the system. 

Although the mean motions n b and n c can be associated with 
the two planets b and c, respectively, it is not the case for the 
secular frequencies g\ and gz, and incidentally, both periastrons 
precess with mean frequency g\ that is retrograde, with a period 
of 474.115 years. The two periastrons are thus locked in an an- 
tipodal state, and the difference Aa> = u> b - to c is in libration 
around 180° with an amplitude of about 36.4° (Fig.©). As a re- 
sult, the argument^ = 2A b - 3A C + lo c librates around 180° with 
the same libration frequency lg. 

4.2. Stability analysis 

To analyze the stability of the nominal solution and confirm the 
presence of the 3:2 reso nance, we performed a global frequency 
analysis (Laska3 fl993h in the vicinity of the nominal solution 
(Fig.© , in the same way a s achieved for the HD 202206 sys- 
tem bv lCorreia et al.l d2005l) . For each planet, the system is inte- 
grated on a regular 2D mesh of initial conditions, with varying 
semi-major axis and eccentricity, while the other parameters are 
retained at their nominal values. The solution is integrated over 
10 kyr for each initial condition and a stability indicator is com- 
puted to be the variation in the measured mean motion over the 
two consecutive 5 kyr intervals of time. For regular motion, there 
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Fig. 5. Evolution of the HD 45364 planetary system over 100 kyr in the co-rotating frame of the inner planet (left) and outer planet 
(right), x and y are spatial coordinates in a frame centered on the star and rotating with the frequency rif, (left) or n c (right). Due to the 
3:2 mean motion resonance the planets never get too close, the minimal distance of approach being 0.371 AU, while the maximum 
distance can reach 1.81 1 AU. We also observe that the trajectories are repeated every 3 orbits of the inner planet and every 2 orbits 
of the outer planet. 



Table 4. Quasi-periodic decomposition of the resonant angle 
6b - 2At - 3 A c + a>i, for an integration over 100 kyr of the orbital 
solution in Table|2] 





Combination 




Vi 


A, 




n h 


n c 


gl 


82 


k 


(deg/yr) 


(deg) 


(deg) 














1 


19.8207 


68.444 


-144.426 








-1 


1 





0.8698 


13.400 


136.931 








1 


-1 


1 


18.9509 


8.606 


168.643 








-1 


1 


1 


20.6905 


8.094 


82.505 








-2 


2 





1.7396 


2.165 


-176.138 








-2 


2 


1 


21.5603 


0.622 


-50.564 














3 


59.4621 


0.540 


-73.279 


1 


-1 








-1 


172.5756 


0.506 


7.504 


1 


-1 











192.3963 


0.501 


-46.923 








-3 


3 





2.6093 


0.416 


-129.207 








2 


-2 


1 


18.0811 


0.420 


121.712 


1 


-1 








1 


212.2170 


0.416 


78.651 





1 


-1 








384.7926 


0.451 


176.155 


1 


-1 








-2 


152.7549 


0.424 


-118.070 





1 


-1 





-1 


364.9719 


0.341 


50.581 





1 


-1 





1 


404.6133 


0.274 


121.729 








1 


-1 


3 


58.5923 


0.212 


-120.210 


1 


-1 








2 


232.0377 


0.201 


24.225 








-1 


1 


3 


60.3319 


0.211 


153.652 


1 





-1 





-1 


557.3682 


0.182 


-86.342 



We have 6), = Yn=i A cos(v, t + 0,-), where the amplitude and phases A ; , 
4>i are given in degree, and the frequencies v, in degree/year. We only 
give the first 20 terms, ordered by decreasing amplitude. All terms are 
identified as integer combinations of the fundamental frequencies given 
in Tablef3] The fact that we are able to express all the main frequencies 
of 9b in terms of exact combinations of the fundamental frequencies gi, 
gi and l e is a signature of a very regular motion. 



is no significant variation in the mean motion along the trajec- 
tory, while it can vary significantly for chaotic trajectories. The 
result is reported in color in Fig. [4] where "red" represent the 
strongly chaotic trajectories, and "dark blue" the extremely sta- 
ble ones. In both plots (Figs.|4^,b), it appears that the only stable 
zone that exists in the vicinity of the nominal solution are the 
stable 3:2 resonant zones. 



1 - 




-1 - 



-1 -0.5 0.5 1 

x(AU) 

Fig. 6. Long-term evolution of the HD 45364 planetary system 
over 5 Gyr starting with the orbital solution from Table[2] The 
panel shows a face-on view of the system, x and y are spa- 
tial coordinates in a frame centered on the star. Present orbital 
solutions are traced with solid lines and each dot corresponds 
to the position of the planet every 100 kyr. The semi-major 
axes (in AU) are almost constant (0.676 < ai, < 0.693 and 
0.893 < a c < 0.902), but the eccentricities undergo significant 
variations (0.123 < e b < 0.286 and 0.042 < e c < 0.133). The 
fundamental periods related to the precession of the periastrons 
are 2nlg\ - -474 yr and 2nlg2 = 3259 yr. 



It is quite rema rkable that, in contrast to the findings of 
ICorreia et all (|2005b for the 5:1 resonance in the HD 202206 sys- 
tem, there is perfect coincidence between the stable 3:2 resonant 
islands, and curves of minimal x 2 obtained in comparing with 
the observations. Since these islands are the only stable zones in 
the vicinity, this picture presents a very coherent view of dynam- 
ical analysis and radial velocity measurments, which reinforces 
the confidence that the present system is in a 3:2 resonant state. 
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In Fig. [5] we plot the evolution of the HD 45364 planetary 
system over 100 kyr in the rotating frame of the inner and the 
outer planet, respectively. Due to the 3:2 mean motion resonance 
trapping, the relative positions of the two planets are repeated 
and never become closer than about 0.37 AU, preventing close 
encounters and the consequent destruction of the system. The 
paths of the planets in the rotating frame illustrate the relation- 
ship between the resonance and the frequency of conjunctions 
with the internal or external planet. The inner planet is in a 3:2 
resonance with the outer planet, so the orbital configuration of 
the system is repeated every 3 orbits of the inner planet and ev- 
ery 2 orbits of the outer planet. In this particular frame, we are 
also able to see the libration of each planet around its equilibrium 
position. 



according to ICrida et all (120081) . an inner disc must be present. 
This singular planetary system may then provide important con- 
straints on planetary formation and migration scenarios. 

The strong gravitational interactions between the planets 
may also allow us to model their effect more accurately in the 
near future. With the current precision of HARPS, about 1 m/s 
for HD 45364, we expect to detect the signature of planet-planet 
interactions in data in a few decades. The planet-planet interac- 
tions may provide important information about the inclination of 
the orbital planes and allow us to determine the precise masses 
of both planets. 

Acknowledgements. We acknowledge support from the Swiss National 
Research Found (FNRS), the Geneva University, Fundacao Calouste Gulbenkian 
(Portugal), and French CNRS. 



4.3. Orbital evolution 

From the previous stability analysis, it is clear that the HD 45364 
planetary system listed in Table|2]is trapped in a 3:2 mean mo- 
tion resonance and stable over a Gyr timescale. Nevertheless, we 
tested directly this by performing a numerical integration of the 
orbits over 5 Gyr using the symplectic integrator SABAC4 of 
lLaskar & Robutell ([2001) with a step size of 0.02 years. The re- 
sults displayed in Fig. [6] show that the orbits indeed evolve in a 
regular way, and remain stable throughout the simulation, which 
corresponds to the estimated age of the star. 

Because of the strong gravitational interactions between the 
two planets, both orbital eccentricities present significant varia- 
tions. The eccentricity of the inner planet is within 0.12 < e\, < 
0.29, while that of the outer planet is within 0.04 < e c < 0.13. 
We also observe rapid secular variations in the orbital parame- 
ters, mostly driven by the rapid secular frequency gi, of period 
Zx/gi « 474 yr (Table|3). These secular variations in the orbital 
elements occur much more rapidly than in our Solar System, 
which should enable them to be detected directly from observa- 
tions. 
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5. Discussion and conclusion 

We have reported the detection of two planets orbiting the star 
HD 45364, with orbital periods of 228 and 342 days, and mini- 
mum masses of 0. 187 and 0.658 Mj up , respectively. A dynamical 
analysis of the system has further shown a 3:2 mean motion res- 
onance between the two planets, which ensures its stability over 
5 Gyr despite the proximity of the two orbits. This is the first 
time that such an orbital resonant configuration has been ob- 
served for extra-solar planets, although an analogue does exist 
in our Solar System composed by Neptune and Pluto. However, 
while Neptune evolves in an almost circular orbit and is much 
more massive than Pluto (which is the largest member of the as- 
teroid family of the Plutinos), the two planets around HD 45364 
have masses comparable to those of Saturn and Jupiter, and are 
evolving in orbits with moderate eccentricity. 

Dynamically, the system is extremely interesting. In the 
nominal solution, the resonant angle 8/, = 2A/, - 3A C + o>b is in 
libration around 0°, with a libration period of 18.16 years and a 
dominating amplitude of 68.44 degrees. Such an orbital config- 
uration may have been reached through the dissipative process 
of planet migration during the early stages of the system evo- 
lution. However, after a resonant capture, subsequent migration 
produces a significant increase in planetary eccentricities, unless 
a damping mechanism is applied. Since the eccentricities of the 
two planets around HD 45364 are relatively small (Table|2]i, mi- 
gration may cease shortly after capture in resonance occurs, or, 



